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ABSTRACT 

We examine the applicability of the stochastic electron acceleration to two high syn¬ 
chrotron peaked blazars, Mrk 421 and Mrk 501, assuming synchrotron self-Compton 
emission of gamma-rays. Our model considers an emitting region moving at relativistic 
speed, where non-thermal electrons are accelerated and attain a steady-state energy 
spectrum together with the photons they emit. The kinetic equations of the electrons 
and photons are solved numerically, given a stationary wave number spectrum of the 
magnetohydrodynamic (MHD) disturbances, which are responsible for the electron 
acceleration and escape. Our simple formulation appears to reproduce the two well- 
sampled, long-term averaged photon spectra. In order to fit the model to the emission 
component from the radio to the X-ray bands, we need both a steeper wave spec¬ 
tral index than the Kolmogorov spectrum and efficient particle escape. Although the 
model provides a natural explanation for the high-energy cutoff of the electron energy 
distribution, the derived physical parameters raise a problem with an energy budget 
if the MHD waves with the Alfven velocity are assumed to be the acceleration agent. 

Key words: acceleration of particles - BL Lacertae objects: individual (Mrk 421, 
Mrk 501) - radiation mechanisms: non-thermal 


1 INTRODUCTION 

Blazars, one of the classes of active galactic nuclei (AGNs), 
have highly collimated je ts with relativistic speed pointing 
close to our line of sight llUrrv fc Padovanilll995l : iMarschen 
I 2 OO 9 II . C haracteristics o f X-ray variability of TeV emitting 
blazars dKataokal I 2 OO 8 II with time-scale of < 1 day mea¬ 
sured in the observer frame suggest that the emitting region 
of the AGN jets that we see in blazars (so-called “blazar 
region”) has the length scale of its structure extending up 
to ~ 10^® cm in its comoving frame when it is estimated in 
this energy band. Since the electromagnetic radiation from 
that region is observed dominantly by virtue of the relativis¬ 
tic beaming effect, blazars are good candidates to investigate 
the energy conversion process of the AGN jets. Photon spec- 
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tral energy distributions (SEDs ) of blazars extend from the 
radio to the gamma-ray bands jFossati et al.lflQ^ l and are 
dominated by non-thermal emission. This demonstrates that 
non-thermal particles are accelerated in the blazar region. 

Blazars are categorized based on their SED peak fre¬ 
quency of the component that extends from t he radio up 
to th e X-ray bands in the vF,j representation (lAbdo et al.l 
I 2 OIOII . For high synchrotron peaked blazars (HSPs), whose 
peak frequency is in the UV or X-ray band, the emission 
spectrum has been primarily interpreted as the synchrotron 
self-Gompton (SSG) process: the synchrotron emission by 
the non-thermal relativistic electrons (and positrons) be¬ 
low the X-ray band and the inverse Compton (IC) scat¬ 
tering of these electr ons off the synchrotron photons for 
the gamma-ray band IIMaraschi . Ghisel lini fc Celottil 1 19941 : 

ISikora k, Madeis5l200ll: fPermer fc Lottll2012ll. 

In modeling observed SEDs, most of the studies assume 
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either a broken power-law form with arbitrary indices for 
the energy distribution of the ra diating non-thermal elec¬ 
trons (e.g., IXavecchio et al.l l200ll 'l or the energy distribu¬ 
tion calculated by solving the kinetic equation under in¬ 
jection of a power-law spectrum, energy losses, and es¬ 
cape (or adiabatic cooling) near and downstream of a shock 
(e.g., iMastichiadis fc KirkI Il997l : iKusunose. Takahara fc Lil 
l200dj . Here we mention three implications on the shape of 
the electron energy distribution (EED) from such studies: 
(i) When a broken power-law EED is applied in the one- 
zone approximation, the change of the spectral index be- 
low and above the bre a k energy tends to be l arger than 1 
jr ^liaferri et aP l2008l : iTavecchio et al.l I 2 OI 0 I : IZhang et al.l 
l2012h . This differs from a prediction of the simple one-zone 
model in which the break is caused by the competition be¬ 
tween radiation cooling and escape, (ii) It is known that the 
broken power-law form of an EED is not always adequate to 
reproduce observed well-sampled SEDs. Purely phenomeno¬ 
logically, adding an additional break t o the EED (i.e., two 
breaks in total) can improve the fit (iKataoka et al.l 1 19991 : 
lAharonian et all l2009l : lAbdo et al.l l2011allbh ~ The physical 
origin of these breaks (or bending of the EED) is not clear, 
(iii) If a population of non-thermal electrons form an SED 
from the radio band (~ Hz) to a higher fre¬ 

quency (e;g;_j_the_STOchrotron peak), a hard EED is re¬ 
quired. IKataoka et al.l (I 2 OO 0 II reproduced the SED of this 
frequency ran ge with the index of the EED of -1.3 5 for PKS 
2155-304, and lKino, Takahara fc Kusunosel (|2002| ~) -1.4, -1.6, 
and -1.8 for PKS 2155-304, Mrk 421, and Mrk 501, respec¬ 
tively. It has been shown that in principle such hard elec¬ 
tron indices can be attained by the first-order Fermi ac¬ 
celeration in the test particle approximation, although the 
results are sensitive to conditions near the shock (e.g., in¬ 
clination angle of the background magnetic field, the shock 
velocity, and the nature of part icle scattering upstream and 
downstream of the shock) (e.g.. iBed narz fc Ostrowska 1 1999 : 
iNiemiec fc Ostrowskill2003 : iLemoine fc Revenu 20061 '). Note 
that since it is uncertain whether an EED of a single popu¬ 
lation is reflected in an SED from t he radio band, it is usual 
not to intend to fit this range (e.g., Krawczvnski et al ]|2004 
IKusunose fc Takahara|[2009 : iRani et al.ll2013f b 

There is also an implication on the acceleration time- 
scale of electrons: (iv) If the balance between synchrotron 
cooling and acceleration time determines the high-energy 
cutoff energy of an EED, the acceleration time of electrons 
with Lorentz factor 7, normalized by the gyroperiod of 
the same electrons, is given as racc(7cut) = 3e/(crTi37cut), 
where e is the charge of an electron, ctt the Thomson 
cross section, B the magnetic field, and 7cut the elec¬ 
tron Lorentz factor at the cutoff. Adopting B = 0.1 G 
and 7cut = 10®, which are derived for emission re¬ 

gions of some T eV blazars with the conventional one- 
zone SSC model dKin o, Takahara fc Kusunosel liool 'l . one 
finds Tacc ~ 10® ( Inoue fc Takaharal Il996l l If it is as¬ 
sumed that only the first-order Fermi mechanism at a sin¬ 
gle shock works in the blazar region, both such a long 
time to accelerate the cutoff electrons and the hard spec¬ 
tra ((iii)) imply that conditions around the shock would 
presumably be very restricted, when extrapolated from the 
test particle simulations (e.g., Bednarz_^_Ostrws^ 1 199a : 
iNieniiec fc Ostrowskil [2004 : ILemoine fc RevemJ 20061 ') . This 
issue may give an important implication on electron acceler¬ 


ation processes, though the spatial dependence of the mag¬ 
netic field and the IC dominance on electron cooling can 
change the value of racc(7cut) to some extent and that one 
could not evaluate Tacc(7cut) well when only a weak con¬ 
straint on the cutoff energy is obtained (e.g., lAbdo et al.l 
l2011alH . 

Instead of interpreting the implications for the non- 
thermal EED mentioned above in terms of the first-order 
Fermi mechanism, stochastic acceleration (SA), which is de¬ 
scribed as diffusion in momentum space, can be invoked to 
be qualitatively consistent with them. Because the models 
including the SA, energy loss processes, and escape from a 
finite acceleration region can form various shapes of an EED 

llBecker. Le fc Dermeil200^ : IStawarz fc PetrokarJlioO^ L the 

SA has the possibility to account for the above features (i)- 
(iv) simultaneously. In addition, since the first-order Fermi 
mechanism necessitates scattering of particles, in the situa¬ 
tion where it really acts slowly as mentioned in (iv), it may 
be accompanied by the SA effectively. 

The emission models including the SA have been 
applied to recent observations of blazars in many 


papers 

('e.g., iKatarzvnski et al.l l2006l: iMassaro et al. 

2001 


lUshio et al.l 2010l: Weidinger. Riiger & Snanier 

201 c 

; 

Lefa, Aharonian & Riegeil 2011: 

Tramacere. Massaro & Tavlorl l201ll: lYan. Zene &; Zhane 

2 OI 2 I: Cao & Wandl2013|). They often focus on either limited 


portions of observed spectra (e.g., near the cutoffs of the 
synchrotron and the IC components) or relatively hard 
flaring spectra, and in some cases add the SA to power-law 
acceleration models or employ multi-zone models. We 
also made a time-dependent simulation wit h a radially 
struct ured jet model including the SA in lAsano et al.l 
ll2014h . In general, hard flaring spectra are compatible with 
a feature of the SA. A possibility of the SA for a slowly 
variable component is im plied with the obse rvation of Mrk 
421 in the X-ray band bv lushio et al.l (120091 '). 

In the present paper, we focus on the issues on the 
shape of EED, which have been raised mainly from one- 
zone steady-state models which assume a (broken) power- 
law EED. In order to investigate the validity of the SA to 
these issues, we apply the one-zone steady-state model em¬ 
ploying a simple form of the SA and SSC emission to well- 
sampled long-term averaged emission of blazars. The obser¬ 
vations compared with the model are the mult iwavelength 
campai gns on Mrk 501 and 421, reported in lAbdo et al.l 
(l2011allHl . 

The paper is organized as follows. We describe our 
model in Section [2] and the procedure of application to the 
observations in Sectional The results are shown in Section 
31 followed by the discussion in Section [S] Section [5] is a 
summary. 


2 MODEL 

We assume that a finite region where electrons are continu¬ 
ously injected into the SA (hereafter referred to as blob) is 
formed intermittently, propagates on a jet, and then reaches 
a steady state, forming the blazar region. We suppose that 
one of the blobs domi nates an observed overall SED at a 
given time on average ilSikora et al.lll99^ . and then calcu¬ 
late the synchrotron and the SSC emission from the blob by 
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solving the kinetic equations of non-thermal electrons and 
photons inside the blob. 

The blob moves relativistically with Lorentz factor T in 
the laboratory frame at an angle to our line-of-sight. 
The Doppler factor of the blob is approximated as T. The 
shape of the blob is taken as cylindrical, characterized by 
the size longitudinal to the direction of its motion and that 
normal to it. Denoting them as i?|| and R± respectively in 
the blob comoving frame, one can write the volume of the 
blob as 


V = ttRIRii . 


( 1 ) 


The two length scales relate to the formation of the blob 
in the jet. Neither acceleration/deceleration nor adiabatic 
change of the blob is considered. 

For simplicity, we assume that the accelerated electrons 
are always isotropically and homogeneously distributed in 
the blob frame and that the time evolution of their phase- 
space distributi on is described as the pure momentum diffu¬ 
sion equation |Hall fc SturrockI 19671: Blaudford fc Eichlerl 


Il987l : lOstrowski fc Siemieniec-Ozieblol 1993). Taking some 
other effects into consideration, we calculate the steady-state 
EE D by solving the kinetic equ ation in the blob frame given 
by JStawarz fc PetrosiarJiiofM i 


dUeij) 

dt 


A 

dj 


2De 


+ 7rad 


UeiA - De 


cine (7) 

dy 


Cte(7) 

^esc,e 


+ 51 ^( 7 ), 


( 2 ) 


where the accelerated electrons are relativistic enough that 
their momentum is approximated as 7meC, rUe is the elec¬ 
tron mass, c the speed of light, ne(7) the number density of 
the accelerated electrons per unit 7, De the diffusion coef¬ 
ficient of the electrons in 7 space due to the SA, 7rad the 
rate of the energy change due to the synchrotron and SSC 
processes, tesc.e the characteristic time for the electrons to 
escape from the blob, and qinj (7) the electron injection rate 
into the acceleration process per unit volume per unit 7. All 
the quantities in equations © and m below are defined 
in the blob comoving frame. The energy distribution of the 
photons emitted by the non-thermal electrons is isotropic 
and homogeneous in the blob frame. Th e photon kinetic 
equation is given by llLi fc Kusunoselbood ) 

fluph(e) . \ I • / ^ I • / ^ Uph(e) 

—^7— = Usynfe) -I- nssA(e) + nic(e) - - -, ( 3 ) 

C'o tesc,ph 

where e is the photon energy in unit of rrieC? , Uph(e) the 
number density of photons per unit e, the first three terms 
on the right-hand side are the photon production and ab¬ 
sorption rate by the electrons per unit volume per unit e, and 
tesc.ph the photon escape time from the blob. The particular 
form of each term in equations m and © is described in 
the following. 

The diffusion of the electrons in 7 space can be caused 
by the scattering off magnetohydrodynamic (MHD) distur¬ 
bances. Assuming that the disturbances are excited uni¬ 
formly throughout the blob and are a superposition of small- 
amplitude waves with the phase speed (in unit of c) of / 3 w = 
Uw/c, we ca n estimate the diffusion coefficient De for fast 
electrons as (iLongairll 19921 : lOstrowski fc Siemieniec-Ozieblol 


ll99'7l:ISchlickeisejl2002|j 



where v {':$> Uw) is the electron speed (« c at present), rg = 
rglA = 'ymeC^/eB the gyroradius of the electrons, B the 
root mean square of the magnetic field randomly oriented at 
scales larger than kinetic ones, diJefr the disturbed magnetic 
field effectively contributing to the scattering of the elec¬ 
trons with Lorentz factor 7 with SB^ff/B being sufficiently 
less than unity, and we consider that the gyroresonance in¬ 
teraction between the disturbances and the electrons at the 
wave number k ~ rg^. The factor {rg/v)~^{5Be«/B)^ rep¬ 
resents the pitch-angl e diffusion coefficient under such sit¬ 
uation (IWentze]||l974h . In order to accelerate electrons to 
higher energies, a wave spectrum distributed in a wide k 
range is needed. Assuming that the wave spectrum is sta¬ 
tionary with a power-law spectrum SB^. oc k~‘‘, containing 
larger amount of energy at smaller wave number, that is, 
g > 1 , we write SB^g as 

(5) 

where Amax is the longest wavelength of the wave spectrum 
and the ratio of the energy density of the disturbed mag¬ 
netic field to the mean magnetic field at the scale of Amax. 
By definition, is less than unity. From equations @ and 
®, one sees De oc 7 '’. 

The escape term describes spatial transport of the par¬ 
ticles in the blob. For simplicity, we take the average length 
scale the particles propagate until escape, denoted by Resc in 
the blob frame, as min{i?x, ^||}- We adopt the spatial diffu¬ 
sion time tesc,e = tesc,e( 7 ) = I?esc/'^ CK 7 '^”^ ^s the electrou 
escape time, where k = lc/3 = rgc{B/SBes)^/9 is the spatial 
diffusion coefficient along the magneti c field, and I oc 7 ^~'^ 
the mean-free path of the electrons llBlandford fc Eichlerl 
Il987l : [Longaiil[r992l L Since treating the spatial transport as 
diffusion is improper for electrons with I > Rese (, in particu¬ 
lar for higher energy electrons if g < 2 ), we set tesc.e = Reac/c 
for electrons with I > Reae in our calculation. As the photon 
escape time, we adopt the light-crossing time Reac/c. 

The rate of the energy change of the accelerated elec¬ 
trons due to radiation can be written as 7 rad = 7 syn + 7 ic, 
where 7 syn and 71 c are the rate due to the synchrotron 
and the IC emission, respectively. The corresponding photon 
source terms in Equation are hayn and hic, respectively. 
Target photons of the IC scattering are only the photons 
originally emitted through the synchrotron process by the 
non-thermal electrons. Synchrotron self-absorption (SSA), 
which is important for the radio band spectra in the ob¬ 
server frame, is taken into consideration only as a sink term 
of photons, denoted by hssA. Electron-positron pair produc¬ 
tion/annihilation is not important for the parameter values 
we will choose and is also neglected. 

The synchrotron loss rate 7 syn, the correspond¬ 
ing photon production rate nsyn(e): s^nd the rate of 
SSA hssA(c) = —cQssAu -ph(c) are calculated following 
iRvbicki fc LightmanI lll979ll and utilizing a numer i cal ap - 
proximation given by iFinke. Dermer fc Bottcheil ||200^ . 
where qssa is the SSA coefficie nt. IC photon production 
rate hic (e) llLi fc Kusunosel[2000l ~l and the corresponding en¬ 
ergy change rate of the electrons 71 c are calculated with an 
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appro xima te form of the IC scatterinR rate, following Ijoned 
(Il968l l and lBlumenthal fc Gouldl (Il97(]l l. 

Without treating thermal population of electrons, we 
just inject electrons to the acceleration process monoener- 
getically at a constant rate, ?inj(7) oc 5(7 — 7inj), where 5 
is the delta function, and 7inj the injection Lorentz factor 
hxed at 10 in this paper. Denoting the total electron injec¬ 
tion luminosity as iinj, we relate ginj with Linj/yinjmeC^F. 

An observed SED of photons radiated isotropically in 
the blob rest frame is calculated with the average photon 
escape rate enph(e)E/tesc,ph, in the vFv representation, by 


vFi, 


e^meC^nph(e)E 

Audi 

^esc,ph 


( 6 ) 



where is the luminosity distance, computed assuming 
the cosmological parameters of Ho = 71kms~^Mpc~^, 
flmo = 0.27, and Hao = 0.73. Photon energy in the ob¬ 
server frame relates to e as hi/ = /{I + z), where 

h is the Planck constant, and 2 the redshift. Steady-state 
energy spectra of the electrons and photons ar e computed 
numerically by solving equa tions Q and ((2| (IPress et al.l 
Il992l : [^ark fc Petrosianlll996l l. We tested our numerical code 
with some an alytical solutions llStawarz fc PetrosiarJ [200^ : 
lMertschll201ll l and found that they are in good agreement. 


Figure 1. Synchrotron spectra calculated for the different value 
of q, the different degree of escape efficiency 7esc(7), and that of 
synchrotron cooling efficiency r^syn (7) . The black thick solid curve 
is for the same value of q, r/esc and rjsyn as adopted in the left 
panel of Figure|2] The black thin dashed (dotted) curve represents 
the case of higher (lower) riesc by a factor of 2, compared to the 
best-fitted one. The blue thick dashed (dotted) curve represents 
the case of higher (lower) 77syn by a factor of 2, compared to the 
best-fitted one. The red thin solid curve represents the case of 
q = 5/3, which is varied from q = 1.85 adopted for the best- 
fitted one. The plotted obs ervational data are the synchrotron 
component of Mrk 421 from lAbdo et al.l (l201ldl and NED. 


3 FITTING PROCEDURE 

To make comparison between observations and the model, 
we give eight physical parameters: P, B, Reac, V, Linj, /3w, 
q, and where the last one is included in equation (Ell- 

Noting that Amax can not be constrained from photon obser¬ 
vations, we do not consider the value of ( or Amax separately. 

The steady-state EED, except for its normalization, 
is determined by the balance among acceleration, cooling, 
and escape; when considering only a synchrotron compo¬ 
nent, one needs to pay attention to only three quantities 
which consist of some of the above parameters: q, syn¬ 
chrotron cooling efficiency rjayn = Vsynil) = tacc/tsyn, and 
electron escape efficiency rjesc = ?7esc(7) = tacc/tesc,e, where 

face — face ( 7 ) ~ 7 jF)a C X 7 tesc.e — tesc.e (7) ^ 7*^ 5 

and fgyn = 'yj'^ayn oc 7“^ (|Stawarz fc PetrosiarJl200i ). 

From equation ©, in the energy range where the ra¬ 
diation cooling can be neglected (?7syn ^ 1), the accelera¬ 
tion forms a power-law EED of Ue(7) oc 7^”"^ at the steady 
state in the energy range of j^esc <C 1, while in the range 
of ?7eac // 1, the escape gradually softens the EED and 
forms a non-power-law one if g < 2. Such effect of parti¬ 
cle escape has the potential to be a natural explanation of 
the shape of EEDs implied by observations mentioned in 
Section [T] ((i)-(iii)). Since the synchrotron component ex¬ 
tends from the radio to the X-ray bands for HSPs, in the 
one-zone model one needs to choose the parameter region 
making electrons with 7 sufficiently lower than 7syn sat¬ 
isfy r^esc > 1, where 7syn is the Lorentz factor defined by 
7syn = 1. Namely, we consider the case that the electron 
escape effect yields a curved EED. Note that the case of 
g ^ 2 is unsuitable for reproducing the well-sampled, grad¬ 
ually softening synchrotron spectra because from equation 
© a power-law EED with the index of dlnne(7)/dln7 = 
1/2—(9/4 -br;e3c)^^^ is formed in the e nergy range of r?svn ^ 1 
for q — 2 (IStawarz fc PetrosiarJlioO^ ) . where r/esc is constant, 


and because for g > 2 higher-energy electrons have shorter 
mean-free path, in other words, shorter ts.ee and longer fesc.e. 
Hence we consider the range of 1 < g < 2 coupled with one 
of the assumptions accompanying equation 

We show in Figure [T] with black curves the variation 
of the synchrotron spectrum with the degree of the elec¬ 
tron escape efficiency. The black thick solid curve fits the 
observed spectrum of Mrk 421 with g = 1.85, rjeaci'y), and 
7syn(7) adopted in the left panel of Figure [5] and displayed in 
the upper panel of Figure [3] The plotted observ ational data 
are the synchrotron component of Mrk 421 fsee lAbdo et al.l 
ll2011bl) for details). The black thin dashed (dotted) curve 
represents the case of higher (lower) ?7esc(7) by a factor of 2, 
compared to the best-fitted curve. The red thin solid curve in 
the same figure represents the case of g = 5/3 with the value 
of ?7esc(7) for the electrons responsible for the synchrotron 
flux at « 10^"^ Hz fixed. This curve (g = 5/3) has stronger 
7 dependence of the escape efficiency than the best-fitted 
curve (g = 1.85) because rieae oc so that the EED 

for the former makes a narrower spectral peak. As can be 
recognized from the black thin dashed and dotted curves, 
if one makes the escape effect less efficient for the red thin 
solid curve, the peak frequency becomes higher (up to which 
the synchrotron cooling becomes effective), and simultane¬ 
ously the EED and so the SED become harder, and vice 
versa. This means that the value of g needs to be closer to 
2, that is, weaker 7 dependence of r;esc(7) than at least the 
case of g = 5/3 in order to obtain SEDs broad enough to fit 
with the observed ones. The adjustment of both the degree 
and the 7 dependence of the escape efficiency is needed to 
explain observed SEDs within this simple one-zone model. 

The blue thick dashed curve and dotted one in Fig¬ 
ure [T] respectively, represent the case of higher and lower 
?7syn(7) by a factor of 2, compared to the black thick solid 
one. These curves show that although the softening of the 
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Figure 2. Comparison between the model SEDs (solid curve) and the observ ational ones (blue ci rcle). The left panel is for Mrk 421, and 
the right panel is for Mrk 501, where the observational data are taken from lAbdo et alJ Hnilbh and lAbdo et alJ ll2011al ') , respectively. 
Non-simultaneous archival data taken from NED are plotted with the gray squares below Hz, for reference. 


SEDs is caused by escape effect, it is insufficient to reproduce 
the high-energy cutoff, and radiative cooling is therefore re¬ 
quired in this case. 

We adjust the model parameters to fit the model to 
the full observed SED including an SSC component, using a 
parameter set obtained from the fitting of the synchrotron 
component as the initial one and considering only the pa¬ 
rameter region satisfying the following two conditions: (1) 
First, the characteristic synchrotron frequency of the elec¬ 
trons satisfying I = J?esc is sufficiently high not to affect 
the observed SED. We consider only the case that the elec¬ 
trons with Lorentz factor 10 times larger than that corre¬ 
sponding to the maximum frequency of the observed syn¬ 
chrotron component satisfy I < 7?esc. This choice means 
that we have the restriction /3w < »?esc' < 1 at this elec¬ 
tron Lorentz factor because rjesc = {I /by dehnition 
and at least ?7esc > 1 for electrons contributing to the flux 
at the frequencies around and above the synchrotron peak 
(as mentioned above in this section). (If) Second, since at 
least the value of {SBes/B)^ for the electrons treated in 
equation m which are responsible for observed emissions 
needs to be sufficiently low, we take {SBeS/B)^ < 0.2 for 
the electrons with the Lorentz factor 10^7syn, which leads 
to < 0.2(rt,(10^7syn))^~'^ from equation 


4 RESULTS 

We apply the model described in Section [5] to the well- 
sampled long-term averaged SEDs of two HSPs, Mrk 421 
and Mrk 501. The results of the fits to the overall observed 
SEDs are shown in Figure [21 The adopted parameters are 
tabulated in Table [T] The observational data for Mrk 421 
was taken from 2009 January 19 (MJD 54 8 50) to 2009 June 
1 (MJD 54983) (Figure 8 in lAbdo et (1201 Ibl E and for 
Mrk 501 from 2009 March 15 (MJD 54905) t o 2009 Au¬ 
gust 1 (MJD 55044) (Figure 8 and Figure 9 in lAbdo et al.l 
ll2011al lL For the gamma-ra y band of Mr k 501, w e adopt the 
data plotted in Figure 9 of lAbdo et al.l (l2011al l with filled 
circles, which do not include the time interval MJD 54952- 


54982 for the time averaging. The excluded term has sig¬ 
nificant flux variations in the gamma-ray band (e.g., a flux 
up to about five times higher than the average one at TeV 
energy) enough to cause a discrepancy between averaged 
Femii-LAT data and non-flaring VERITAS data in overlap- 
ping energy range if included for the time averaging (see 
lAbdo et al.l (l2011ah for more details). During the campaign 
Mrk 421 was in a low activity state at all wavebands (e.g., 
flux variations are typically less than a factor of two relative 
to the average flux level in gamma-ray band) (lAbdo et al.l 
l2011bf ). The data we adopted are suitable ones currently 
available for modelling the typical low state of these objects. 
The archival data plotted in Figure[2]for reference was taken 
from NED. 

As shown in Figure (2] the two samples appear to be 
fitted well. Here we should note that the reduced chi-square 
values of the fits shown in Figure [2] are 3.9 and 1.3, re¬ 
spectively, for Mrk 421 and Mrk 501 for the gamma-ray 
band (3.8 x 10^^-9.6 x 10^® Hz for Mrk 421 and 3.9 x 10^^- 
1.1 X 10^^ Hz for Mrk 501)13, though they are not minimum 
ones. The current work is not able to specify the causes of 
the statistical poorness of the fits. This may possibly indi¬ 
cate the need for non-trivial spatial structure of the emis¬ 
sion region and/or additional emission components such as 
an external Compton component in order to reproduce the 
actual data more accurately. As an example for the latter, 
the comptonization of an external radio photon field has 
been proposed so as to explain somewhat but systemati- 
cally undere s timat ed flux near 10^® Hz (see Section 4.2 in 
lAsano et al.l ll2014l') j. Nevertheless, we consider that the fits 
are sufficient for investigating the validity of SA to the issues 
described in Section [T] 

The characteristic synchrotron frequency of the elec¬ 
trons with the Lorentz factor 7inj is about 10® Hz for Mrk 
421, and about 10® Hz for Mrk 501, measured in the observer 
frame. The break seen near 10^® Hz in the both panels in 

^ For the frequency range below the X-ray band, since the ob¬ 
served data show different flux levels at very close frequencies with 
small statistical errors, we do not evaluate chi-square values. 
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Table 1. Physical parameters used in the calculation shown in Figure [2] 



r 

B 

Resc 

V 

Tinj 

/3w 

Q 

^max 

7inj 

unit 


G 

lO'^^cm 

10^' cm'^ 

1039erg 



lO-l^cml-? 


Mrk 421 

35 

0.081 

3.5 

1.9 

1.6 

0.16 

1.85 

4.7 

10 

Mrk 501 

31 

0.014 

28 

95 

15 

0.23 

1.92 

0.12 

10 


this figure is caused by SSA. In the upper panel of Figured 
the characteristic times of acceleration, escape, and cooling 
are displayed as a function of 7 for the parameters adopted 
for Mrk 421. The lower panel shows the calculated FED of 
Mrk 421. It can be recognized from this fignre that grad- 
nal softening from the radio band to higher frequency is 
formed by an increasing dominance of escape over accelera¬ 
tion with an increasing electron energy and that the cutoff 
at the synchrotron peak is cansed by the balance between 
acceleration and cooling time under rapid escape. Note that 
not but mimics the synchrotron spectrum 

dTramacere. Massaro &: Tavloijlioill l. Not only the soften¬ 
ing by the efficient escape, but also cooling is required to 
well reproduce the observed shape of the cutoff in this case. 

The synchrotron component of Mrk 501 has a softer and 
broader portion, compared to Mrk 421, as seen in Figure[2l If 
one interprets such an SED from the radio band in the one- 
zone model, q closer to 2 (adopted value is 1.92) is needed 
as mentioned in Section [S] Treating the slope of the wave 
spectrum as a free parameter is essential to apply our simple 
model to various objects. At least the phenomenologies such 
as Kolmogorov-type q = 5/3 and Krai chnan-type g = 3/2 
dColdstein. Roberts fc Matthaeu^ll995ll are not fit for the 
application here; larger q but smaller than 2 is needed. We 
omit the figures for the time-scales and EED for Mrk 501 
because it is basically similar to that for Mrk 421. 

The adopted parameter values are listed in Table [T] We 
call them and the ones obtained with them the canonical 
values. By examining the relation between a model SED 
and physical quantities, here we show that some of them 
are constrained to the canonical values with the SED fitting 
while others are not. First, we note that the synchrotron 
component is not modified except by the influence of the 
IC cooling when we change the parameter values with the 
flux normalization of the synchrotron component and the 
values of r^esc and ?7syn fixed as a function of synchrotron 
frequency measured in the observer frame, i^syn, with the 
corresponding electron Lorentz factor 7, i.e., we obtain 


* 

^syn 

= i = (B7^r)*, 

(7) 

1 \ * 
^'/syn 

= 1 = (TsynT^)* = {B^-yUaV^Vy , 

(8) 

Vsyn 

= i = (bV^/3-=(Ca)/;J)-1)‘, 

(9) 

hyn) 

= 1 = V')% 

(10) 


where the ratio of a quantity to the canonical value is de¬ 
noted by the superscript *, (!zFh)*yn is the ratio of the 
flux normalization of the synchrotron component, Lgyn the 
synchrotron luminosity measured in the blob frame, ub = 
/8n the energy density of the mean magnetic field, Ue that 
of the non-thermal electrons, and L*y^ = (ubJUbV)*. Next, 
for the fixed synchrotron component, the IC component is 


Table 2. Energy contents in the blob for the parameters shown in 
Table[T] the energy density of magnetic field, ub \ that of photons, 
iipli; and that of non-thermal electrons, u^. 


Ub Mph Uf, 

unit 10“^ erg cm“® 10“^ erg cm“® 10“^ erg cm~® 

Mrk 421 0.03 0.02 9 

Mrk 501 0.0007 0.002 2 


not modified when: 


(Usyn / Ub ) 

7* 

F* 


1 = (yUefllesc)* , 

1 , 

1 , 


( 11 ) 

( 12 ) 

(13) 


where Msyn is the energy density of the synchrotron photons, 
measured in the blob frame, and 1/*^^ = {usynV/Reac)*■ Note 
that the frequency range of the Thomson scattering regime 
of the IC component relative to the synchrotron component 
and the IC cutoff frequency roughly scale as 7*^ and (7r)*, 
respectively. From equations ©-(Ha}, we obtain 



= 7*=r‘ = i, 

(14) 

* 

lie 

= 

(15) 

/?* 

-^esc 

= 

(16) 

/?: 

= 

(17) 

l-qy 

'max/ 

= V*~^. 

(18) 


Thus, B, F, and the total electron energy UeV are well con¬ 
strained to the canonical values. Moreover, the SED is un¬ 
changed if the physical parameters of the blob are changed 
according to equations (P ^ -(fT5 |) with V*. That is, a param¬ 
eter set adopted for an observed SED as shown in Figure [2] 
includes only one degree of freedom and can vary the pa¬ 
rameter values within the restrictions (I) and (II) in Sec¬ 
tion [3] Although Linj, one of the fitting parameters, does 
not appear in the above analysis explicitly, one can rewrite 
equation dm in terms of Tinj. Note that all the presented 
figures are independent of the uncertainty and our general 
conclusions are unaffected. 

For concreteness, we have shown in Table[T]the parame¬ 
ter values of Resc, V, Tinj, and by adopting the maxi¬ 

mum allowed value of /3w from the restriction (I) described in 
Section|3] The shown parameter sets can be scaled by V* in 
the range of 0.04 ^ V* ^ 1 for Mrk 421 and 0.09 < "F* sS 1 
for Mrk 501 with the SEDs fixed, where the upper and lower 
bound come from (I) and (II) in Section [3] respectively. 

Energy densities of respective components, measured 
in the blob frame, are listed in Tabled Although the mag¬ 
nitude relation between a nd mb is consistent with the 
conventional SSC model le.g.. [lGno. Takahara fc Kusuno^ 
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Figure 3. Upper panel: The characteristic times of acceleration, 
escape and cooling for non-thermal electrons, corresponding to 
the left panel of Figure [2] The red thick solid line represents the 
acceleration time face, the green dashed line the electron escape 
time tesc,e, the blue thin solid curve the radiation cooling time 
at the steady state, 7 / | 7 rad|! fhe cyan dotted line the syn¬ 
chrotron cooling time tsyn- Lower panel: Electron energy distribu¬ 
tion of the non-thermal electrons, multiplied by 7 ^. The electrons 
are injected at 7 = 10 continuously at a constant rate. 

I2002I1 . the non-thermal electrons possess much greater en¬ 
ergy than the magnetic field, rte/ws ^ 1, compared to the 
previous results. This difference originates from relatively 
small f?esc, that is, the importance of the escape effect as 
seen from equations m and ca. Within the uncertainty 
of the parameters, we can only make Ue (and Ub/ub) larger 
from the canonical values since (ueV)* = 1 is required to fix 
the SEDs from equation (1151) . where E* ^ 1. 

5 DISCUSSION 

5.1 Shape and the Number of the Blobs 

As seen in Table [T] a relatively small Reaci<^ is 

required in order to make the escape effect work to form 
the gradually softening SEDs. This means that the emission 
region is either disk-like or spindle like. 

However, the extreme shape of the blob may be alle¬ 
viated if one assumes that a number of the blobs with the 
same properties simultaneously contribute to an observed 
spectrum. The number of the blobs can be easily in¬ 
cluded into the analysis in Section |4] by multiplying Layn in 
equation dll) by Nh- Then V* in equations GSl) - dll is re¬ 
placed with (VNh)*■ It can be recognized from equation fT]) 
and the equation modified from equation GSl) with this re¬ 
placement, Alesc = (VNh)*, that for example, when increas- 



hcghini) 

Figure 4. The time-scale of damping of MHD disturbances 
through acceleration of non-thermal electrons in unit of l/fcuw 
for the case of u^{k) ~ SB^f^/Sir. The horizontal axis represents 
wave number k normalized by the inverse of the gyroradius of the 
electrons with the injection Lorentz factor 7inj. The dotted line 
indicates a lower limit of the wave cascading time. 

ing A7 and decreasing V with (VNb)* = 1, one can scale 
the ratio iiesc/max{i?_L, i?||} = min{i?x, i?||}/max{i?x, J?|| } 
up to unity as if R± > = Reac or 

as = Nh if J?|| > R± = Resc without changing 

both the SED and the parameter values (except for Linj). 
This is the case of multiple blobs contributing to an SED. 
In the above example, we have i?esc/inax{i?x, I?||} < 1 with 
Vb < 10® and < 10^ for Mrk 421 and Mrk 501, respectively. 
However, a large number which leads to R± ~ i?|| is rather 
unlikely in the spirit of the simple one-zone model approach. 

5.2 Energy Budget 

Let us consider the energy balance between the waves and 
the accelerated particles. Since we have supposed the net 
energy transfer between the waves and the non-thermal elec¬ 
trons is from the former to the latter, the energy supply from 
large scale disturbances to small scale ones has to be rapid 
enough to meet the assumption of the stationary wave spec¬ 
trum. We estimate the time-scale of damping of the wave 
energy caused by the electron acceleration, as a function of 
k, as 

idau.(fc) ~ ^ , (19) 

L/acc 

measured in the blob frame, where u^{k) is the to¬ 
tal energy density of wave modes around k, Wacc = 
7®meC®ne(7)/tacc(7), and k ~ r“®. If this is too short, 
waves would damp through particle acceleration and then 
Mw(fc) would be modified. For illustrative purpose we dis¬ 
play tdam(fc) in Figure [7] for the steady state we have ob¬ 
tained, supposing Ww(fc) ~ 6Baff/8n. By comparing tdam(fc) 
to 1/fcuw, which represents a lower limit of the time-scale of 
the wave cascading, we find that tdam is smaller than l/kv„ 
in an intermediate range of k. This means that before the 
blob reaching the steady state, the assumed wave spectrum 
would be modified and evolve toward a different equilibrium 
state. The model is hence inconsistent in this case. 

Under our treatment of the electron transport, this en- 
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courages us to interpret Uw as the sound speed (not as the 
Alfven velocity) since the wave energy can be much larger 
than the magnetic part in this case if the plasma beta of 
the emitting region is much larger than 1. It seems qualita¬ 
tively reasonable because the energy of the magnetic field 
is much less than that of the particles in the blob. This 
implies that the fast mode waves may be the acceleration 
agent. (In this case u^{k) and so tdain{k) become larger in 
proportion to the plasma beta although the cascading time 
is also affected.) The parameter sets shown in Table [T] have 
the value of Uw close to the upper limit of the sound speed 
(i.e., c/y/Z), which seems likely in the blazar region. Any¬ 
how, more sophisticated turbulence modeling accounting for 
particle acceleration is needed for self-consistency, which is 
beyond the scope of this paper. 


6 SUMMARY 

Motivated by the previous implications about the energy 
distribution of the non-thermal electrons in the blazar re¬ 
gion ((i)-(iv) described in Section [T]), we have examined the 
stochastic acceleration by a weak random electromagnetic 
field as the acceleration mechanism of the electrons by mod¬ 
elling the SEDs of two HSPs, Mrk 421 and Mrk 501, from the 
radio to the gamma-ray bands. All the model parameters are 
time independent (Section [2|. The long-term time averaged 
SEDs appear to be well fitted by the emitting blob reaching 
a steady state (Figure [21) • The blob can be considered to 
represent the jet emission region(s) mainly contributing to 
the observed SED on average in the long term. The good¬ 
ness of the fits in gamma-ray band is, however, statistically 
not sufficient as mentioned in Section [4] and more refined 
modelling are needed to reproduce the actual data more ac¬ 
curately. With the current simple model, we have obtained 
the following results. 

The model can naturally produce the broad and curved 
electron energy distribution, which have been suggested 
based on the conventional models so far (Section [T|). To re¬ 
alize such electron energy distribution, two points are cru¬ 
cial in the model: a stationary, power-law spectrum of MHD 
waves with a steeper spectral index than the Kolmogorov 
spectrum at the relevant wavelengths and efficient particle 
escape from the acceleration region (Figures [T] and O). The 
former is associated with the latter through weaker energy 
dependence of the escape efficiency. As for the formation of 
the high-energy cutoff of the synchrotron component, radia¬ 
tive cooling becomes important. 

These results lead to a few issues with the obtained 
parameter values. First, the shape of the blob is very asym¬ 
metric having a relatively small side length compared with 
the volume of the blob to realize an efficient escape. Sec¬ 
ond, since emitted photons also efficiently escape, in order 
to have a reasonable value for the ratio of the energy den¬ 
sity of the synchrotron photons to that of the magnetic field, 
which is strictly constrained by an observed spectrum, our 
model requires a large ratio of the energy density of the 
non-thermal electrons to that of the magnetic field (equa¬ 
tion GU), which ratio is actually orders of magnitude larger 
than the ones hitherto adopted for the two HSPs (Section 

m. 

A trivial but important issue of the model is whether 


and how such properties phenomenologically required for the 
emission region can be formed in the blazar jets. If we esti¬ 
mate the energy density of the MHD waves responsible for 
the electron acceleration on the assumption that they prop¬ 
agate at the Alfven velocity, the derived physical parameters 
raise a problem with an energy budget of MHD turbulence, 
i.e., the turbulence would be damped effectively by the par¬ 
ticle acceleration fSection l5.2ll . The model is inconsistent in 
this case. The difficulty may be reduced to some extent if 
we assume that the waves propagate at the speed of sound 
in a high-beta plasma (i.e., the fa st mode waves). 

Recently, lAsano et al.l (l2014l l considered the emission 
from the electrons stochastically accelerated in the contin¬ 
uous and radially-st ructured jet , wher e the escape effect 
can be negligible. In lAsano et al.l (l2014fl assuming the Kol¬ 
mogorov spectrum [q = 5/3) for the turbulent magnetic 
field, which is different from our treatment of q, they in¬ 
vestigated the effect of the radial evolution of the physical 
parameters on the electron energy distribution and observed 
SED. They have shown that electron spectra become soft as 
required even when q = 5/3 by assuming the evolution of 
the electron injection rate, that of the momentum diffusion 
coefficient, and adiabatic cooling. Although their jet model 
is quite a contrast to the blob model in this paper, it is also 
able to form the curved electron energy distribution that can 
explain the same SED of Mrk 421 by the stochastic acceler¬ 
ation. Thus, although the stochastic acceleration is a viable 
model, the way to realize a broad but a non-power law elec¬ 
tron spectrum is not unique and further work is needed to 
solve the problem. 
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